Spin-chirality duality in a spin ladder with four-spin cyclic exchange 
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Effect of four-spin cyclic exchange on magnetism is studied in the two-leg S — 1/2 ladder. We 
develop an exact spin-chirality duality transformation, under which the system is self-dual when the 
four-spin exchange J4 is half of the two-spin exchange. Using the density-matrix renormalization- 
group method and the duality relation, we find that the four-spin exchange makes the vector chirality 
correlation dominant. A "chirality short-range resonating-valence-bond" phase is identified for the 
first time at large J 4. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg, 74.25. Ha 



Recently, it has been realized that the two- leg S — 1/2 
ladder compound La :E Cai4_ a ;Cu24O4i0,|^0j4|] and two- 
dimensional (2D) antiferromagnet I^CuC^UQ have a 
certain strength of four-spin cyclic exchange interactions. 
Theoretically, four-spin cyclic exchange emerges in the 
strong-coupling expansion of the one-band HubbardQ 
and d-p L S] models as the leading correction to the nearest- 
neighbor two-spin exchange. Cyclic exchanges were also 
found to be large in magnetism of 2D quantum solids, 
e.g. solid 3 He filmsQ and Wigner crvstalsfioj. The effect 
of four-spin cyclic exchange on magnetism is, however, 
hardly understood, since it has frustration by itself: the 
question of what type of magnetic ordering tends to be 
realized by the four-spin exchange is still unsettled. For 
example, in the context of magnetism of solid 3 He films, it 
was argued that the four-spin exc hang e on the triangular 
lattice can induce scalar chirality though finite-size 
system analysis could not find evidence for such ordering, 
instead showing spin-liquid ground states[l2j|. 

To clarify magnetism induced by the four-spin cyclic 
exchange, we consider the spin ladder. Spin ladder anti- 
ferromagnets have been attracting extensive interest be- 
cause they have a spin gap, a short-range resonating- 
valence-bond fRVB) ground state, and superconductivity 
upon doping[l3j. On the two- leg 5 =1/2 ladder it was 
shown numerically that the spin gap decreases rapidly 
with increasing the four-spin cyclic exchange J4Q and 
vanishes at a critical point, J4/J ~ 0.3 where J is the 
two-spin exchange [l4l Ho| . The nature of the new phase 
for large J4 was not established. 

In this letter, we show that the four-spin exchange in 
the two-leg 5=1/2 ladder has a tendency to induce 
a vector chirality correlation. First, we describe an ex- 
act duality transformation between the Neel-spin opera- 
tor (s u — s 2 .f) /2 and the vector-chirality one x s 2 j 
on the rungs, where s Mi / is the spin operator at the 
site on the leg fi = 1,2 and the rung I, The system 
is self-dual under this transformation at J4/J = 1/2, 
where the Neel-spin and the chirality interchange their 
roles: the former gives the most dominant correlation 
for small J4 while the latter does for large J4. Us- 



ing the density-matrix renormalization-group (DMRG) 
method 16] and the spin-chirality duality, we studied the 
ground-state phase diagram of the ladder for the whole 
region of < J4 < 00. We find the spin short-range 
RVB phase, an intermediate phase with a very small spin 
gap, and a novel chirality short-range RVB phase. Our 
findings of exotic magnetic states with dominant vector- 
chirality correlation at large J4 suggest that the four- 
spin exchange can induce exotic electronic states in doped 
systems such as high-T c superconductors, whereas only 
two-spin exchanges have been taken into account in t-J 
models in searching the mechanism. 

Let us consider the Hamiltonian defined as 



n 



1 



{JmngSl,; ' S 2 .l + Jleg(Sl,i ■ Si j+l 



"S2,i ' S 2 ,i + l) + J4K1}, 



(1) 



where J rU ng (J\eg) denotes the two-spin exchange con- 
stant on rungs (legs) and Ki the four-spin cyclic exchange 
on a plaquette {(1, 1), (2, 1), (2, 1 + 1), (1, 1 + 1)}, 



Ki = Si,j • 82,1 + Si,; + i • s 2 ,i+i + 

+ S2 ) cS 2l l + l+Si,i'S2,J+l + 

+ 4(si,j • s 2 ,;)(si,;+i ■ s 2,i+i) 

+ 4(si,( • Si,j4-l)(S2,J • S 2 ,i+l) 

- 4(si,j • s 2ji+ i)(s 2ji ■ sij + i). 



Sl,/ • Sij + i 
$2,1 ■ Sl,i+1 



(2) 



All the coupling constants are assumed to be positive, 
>Aung, J\cg, J a > 0. It is instructive to rewrite the Hamil- 
tonian (JJJ as 



H 



(Jr, 



•2J 4 ) ^Sl,! 



s 2,l 



(Jleg + J4) S Mi ; • S Mi ; + i 

J4 (Su • S 2: ; + i + Sij + i • B2,j) 
I 

4J4 (si,j • s 2 ,i) (si,;+i • s 2 ,;+i) 
1 

4J 4 ^ ( S U X S 2,0 • ( S l,2 + 1 X S 2,i + l) ■ 



(3) 
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An interesting contribution of the cyclic exchange ap- 
pears in the last term, which introduces a coupling be- 
tween vector chiralities on nearest-neighbor rungs. This 
term tends to induce non-zero vector chiralities on every 
rung arranged in an antiparallel pattern. Hence, one can 
naively expect that for large J4 the vector chirality be- 
comes an important degree of freedom although the frus- 
tration between the various terms in eq. © complicates 
the situation. We will show later that the vector-chirality 
correlation indeed becomes dominant for large J4. 

To elucidate the relation between the spin and chirality 
degrees of freedom, we construct a duality transformation 
between them. Let us begin with the commutation rela- 
tions between the total rung-spins W; = si,; + s 2 j and 
the vector chirality V; = 2 Si,; x s 2 j given by 



v{*,v£ 



0/371/7 ; 



We note that the commutation relations are identical to 
those which hold between the angular momentum and the 
Runge-Lenz vector of an electron system in a hydrogen 
atom. We can disentangle the algebra by introducing 
new operators defined by 

S i = \ ( W * - v = \ ( 8 i.l + s 2-0 - s hi x s 2,h (4) 

T ( = ~ (W, + V;) = i (81,1 + S 2 ,l) + 8i,j X 82,1. (5) 

These operators obey the commutation relations, 



If, if, 



= 0, 



and satisfy = T? = 3/4. Thus, the new opera- 
tors S/ and T; are S = 1/2 pseudo-spin operators de- 
coupling each other. It is interesting to note that the 
original spins Si,; and s 2 ,; may be expressed in terms 
of Si and T/ simply by interchanging their roles in 
eqs. gj| and ©, i.e., Si,j = § (S/ +T/) + S ( x T; and 
S2,i = I (S; + Tj) — Si x T ; . We hence call this a "dual- 
ity" transformation. The relations between the original 
and new spin operators are summarized as 



Sl,i + s 2 ,z 
si,i - s 2 ,; 
-2 si,i x s 2 .i 
si,i • s 2 ,; 



2 Si x T h 

S;-T ; , 

S/-T ; . 



The transformation thereby exchanges the Neel-spin and 
the vector chirality on the same rung. 



In terms of the new spin operators, the Hamiltonian 
(JSJ is rewritten as 

H = ( Jrung + 2 J 4 ) ^ Si ■ T; 

I 

^£ + 2J^53(S l -S I+ i+T I -T I+1 ) 
+ ^i^( Si .T i+1 +T ; .S i+1 ) 

! 

+ 4J 4 ^(S r T ; )(S w -T m ) 

i 

+ 2J lcg (Si x T ' (Sj+i x T i+ i) . (6) 
1 

Thus, the duality transformation leaves the form of the 
Hamiltonian unchanged and only affects the coefficients 
of the second, third, and fifth terms. An interesting ob- 
servation here is that in the case of J4 = J\ cg /2 the origi- 
nal Hamiltonian Ti and its dual TL are equivalent includ- 
ing the coefficients. Hence the Neel-spin (si,; — S2,;) /2 
and the vector chirality si,; x s 2 .z show identical behavior 
on this "self-dual" line. 

To clarify the consequence of the spin-chirality dual- 
ity around J4 = J\ cg /2 and the nature of ground states, 
we study the low-energy properties of the Hamiltonian 
@ numerically. For simplicity, we focus on the case of 
>Aung = ^lcg = J and investigate the ground-state phase 
diagram on the J 4/ J line hereafter. Using the DMRG 
method^], we have calculated the energy gap of spin 
excitations 

A 0M (L) = E (L; M) - E a (L: 0), (M = 1, 2), (7) 

where Eq(L;M) is the lowest energy in the subspace of 
s totai = 12i( s i 1 + s 2 1) = M m a, finite ladder of L rungs. 
For the best performance of the DMRG method, an open 
boundary condition was imposed. We have also calcu- 
lated the ground-state spin correlation functions 

Co(r) = j<K* (8) 

c'M = J<«*-40«*'-4*')>* (9) 

and the vector-chirality correlation function 

C K (r) = (( Bl>1 x s 2 ,;) z (si,,/ x sa,,,)*), (10) 

with I = l - r/2 and I' = l + r/2 17J. The index l 
represents the center position of the open ladder, i.e., 
l a = L/2 for even r and Iq = (L + l)/2 for odd r. We 
have employed the finite-system method with improved 
algorithm[l|j| and kept up to m — 500 states per block. 
The numerical errors due to the truncation are estimated 
from the difference between the data of different m's. The 
system size is up to 2 x 100 sites. 
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FIG. 1: System-size dependence of the spin gaps for (a) 
(J, Ja) = (1,0.1); (b) (1,0.3); (c) (1,0.5); and (d) (0,1). The 
gaps Aoi(-L) and Ao2(L) are plotted by circles and squares, 
respectively. The numerical errors of the DMRG calculation 
are smaller than the symbols. 



First, we discuss the parameter region J4/J < 1. It 
is known that for J4 = the system belongs to a spin- 
liquid phase, in which the ground state is well described 
by the RVB state. The spin gap is open in this phase. 
We show in Fig. ^our numerical results for the spin gaps 
Aqm{L). The data are extrapolated by fitting them to a 
polynomial form, Aqm(L) = A m(oo) +a/L + b/L 2 . For 
J4/ J < 0.3, the spin gaps decrease smoothly for both M's 
as L increases, and consequently, the extrapolation works 
pretty well [see Fig.^a) and (b)]. The extrapolated val- 
ues Aqm(oo) are shown in Fig.[2with J = 1 in the region 
J 4/ J < 1. The spin gaps decrease smoothly as J 4 in- 
creases from and vanish around J4/J ~ 0.3, suggesting 
a phase transition accompanied by vanishing of the spin 
gap at J4/J = {Ji/J)d — 0.3. Unfortunately, accurate 
estimation of the critical value (J4/ J) c i is quite difficult 
due to the very slow vanishing of the spin gaps around the 
transition point. When J4/J is larger than 0.3, Aoi(L) 
shows bumpy behavior as seen in Fig. ^c). This may 
be attributed to effects of open boundaries. The value 
of Aoi(L) becomes exactly within numerical accuracy 
for several L, which suggests a spin-triplet ground state. 
On the other hand, the spin gap Aq 2 (L) exhibits a rather 
smooth L-dependence even for J 4 / J > 0.3. The extrap- 
olated gaps Ao2(oo) are very small, less than 0.02, but 
seem to be finite for this region. Very recently, Lauchli 
et al. studied independently the same model but with 
a different boundary condition and showed that the sys- 
tem for this parameter region belongs to a phase with 
a very small gap exhibiting the translational symmetry 
breaking 0. The work of Ref. also pointed toward 
this result. Our results are thus consistent with theirs 
although the number and type of excitations in the finite 
systems might differ from each other because of the dif- 
ferent boundary conditions. Further studies, especially 
by analytic methods, are desirable for clarifying the na- 
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FIG. 2: Extrapolated spin gaps in the limit L — > 00 as func- 
tions of J4 (left), J (center), and J IT (right); see text. The 
circles and squares represent Aoi(oo) and Ao2(oo), respec- 
tively. Inset: Enlarged figure for 0.3 < J4/J < 2. The error 
bars represent those from the extrapolation procedure. 



ture of excitations and why the spin gap is so small in 
the entire region of the phase. 

Figure [3] shows the spin correlations Cg(r) and C*(r), 
and the vector-chirality correlation C K (r) for several typi- 
cal sets of parameters. For J4/J < (J4/ J) c i — 0.3, all the 
correlations decay exponentially, which is consistent with 
a finite spin gap. The Neel-spin correlation C£(r) is the 
strongest among the calculated correlation functions [Fig. 
OJa)], as in the usual antiferromagnetic (AF) ladder. For 
J4/J > ( J4/ J) c i, on the other hand, the correlation func- 
tions decay very slowly reflecting the small spin gat»|2o|. 
We find that the Neel-spin correlation C*(r), which is 
dominant for small J4, keeps reducing as J4 increases 
while the vector-chirality correlation C K (r) keeps grow- 
ing and becomes dominant for large J4. They interchange 
with each other at the self-dual point J4/J = 0.5; we can 
see their perfect coincidence in Fig. E^c) . These results 
strongly suggest that the system exhibits symmetric be- 
haviors with respect to the self-dual point J4/J = 0.5 
with exchanging the roles of the Neel-spin and the vector- 
chirality correlations. We also note that the dimer opera- 
tor sxj-sij+i— S2J-S2J+1 and the scalar-chirality operator 
(si,; + «2,0 • x £2,1+1) + (si.i x s 2t i) ■ (si,;+i + s 2 ,i+i) 

are related to each other by the duality transformation, 
and consequently, their correlation functions must inter- 
change exactly at J4/J = 0.5, which is consistent with 
numerical results in Ref. Il8t 

Next, we consider the parameter region J4/J > 1. 
Hereafter, we set J4 = 1. It can be seen in Fig. [21 that 
the spin gaps open for J4/J > (Ji/J) C 2 — 1/0.4 = 2.5. 
Again, accurate estimation of (Ja/J)c2 is quite difficult 
due to the very slow opening of the gaps. We also note 
that for large J4/J the spin gap A01 exhibits a smooth 
L-dependence and does not become for finite L [see Fig. 
ID!-.)], suggesting the absence of the triplet ground state in 
the finite systems. To elucidate the nature of the system 
in this large J4/J region, we consider the case J4/J = 00 
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(a) J = 1, J 4 = 0.1 » (c)J = 1,J 4 =0.5 




(b)J = l,J 4 =0.3 » (d) J-1, J«-1.0 




FIG. 3: Correlation functions (r) (circles), |Co(r)| 
(squares), and |C K (r)| (triangles) as functions of distance for 
J = 1 and (a) J 4 = 0.1; (b) J 4 = 0.3; (c) J 4 = 0.5; and (d) 
J a = 1. The system size is L — 80. The error bars represent 
the numerical errors of the DMRG calculation. 

(J4 = 1 and J = 0) using the spin-chirality duality trans- 
formation. In this case, the transformed Hamiltonian © 
is expressed as 

7~Lj4 = Xung ^ S; • T; + Ji cg ^ (S; ■ S; + l + Tj • T; + l) 

I I 
+ J rr ^(S i -T ; )(S i+1 -T i+1 ) (11) 

with J rung = ./log = 2 and J rr = 4. Notice here that, if 
one sets J rr = in eq. fTTJ) . the system is equivalent to the 
usual two-leg AF ladder, which has the short-range RVB 
ground state consisting of the spins S; and T;. In Fig. 
we show the J rr -dependence of the extrapolated spin gaps 
A m(oo) for < J lr < 4. It is clear that the spin gaps 
remain finite for the entire region of < J rr < 4 and are 
smoothly connected to the spin gaps at J rr = 4; there is 
no phase transition between J rr = and J rr = 4. We thus 
conclude that the Hamiltonian l|ll(l with J rr = 4, and 
accordingly, the original ladder Q in the limit J4/ J = 00 
belong to the same RVB phase as the AF ladder of the 
spins Si and T; with J rr = 0. Small size of the spin 
gap at J4/J — 00 can be understood from the fact that 
the system is close to the quantum critical point between 
the short-range RVB and intermediate phases. Note that 
the dominant correlation function in this RVB phase is 
that of the Neel-spin (S/ — T;) /2 and hence, in terms of 
the original spins, the correlation of the vector chirality 
Si,i x s 2 ,i is the strongest. We therefore term this novel 
phase the chirality short-range RVB phase. 

To summarize, using the spin-chirality duality trans- 
formation, which is developed in this letter, as well as the 
DMRG method, we have found that the four-spin cyclic 
exchange makes the vector chirality correlation domi- 



nant. The chirality RVB phase appears for large J4. It 
has been found that the system exhibits symmetric be- 
havior with respect to the self-dual point J4/J =1/2 
by interchanging the Neel spin and the vector chiral- 
ity. We remark that the duality transformation is ap- 
plicable to any spin Hamiltonian, since it is based only 
on the spin commutation relation. This transformation 
should be useful in studying various topics. One exam- 
ple is the spin-orbital model around the SU(4) symmet- 
ric point [2l|. We have found in the two-leg ladder with 
extended four-spin exchange that the self-dual line con- 
nects with the SU(4)-symmetric point Another ex- 
ample is a magnetization plateau induced by the four- 
spin exchange[23j. Since the total spin i is in- 
variant under the dual transformation, the duality rela- 
tion holds even in a magnetic field. Effect of four-spin 
exchange on hole-doped systems is also to be considered. 
It would be interesting to investigate relation to possible 
hidden orders proposed for high-T c superconductors, e.g. 
the staggered currents. 
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